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Abstract 

We consider the production of the Standard Model Higgs boson through the gluon 
fusion mechanism in hadron collisions. We present the next-to-next-to-leading 
order (NNLO) QCD result of the hard-collinear coefficient function for the all- 
order resummation of logarithmically-enhanced contributions at small transverse 
momentum. The coefficient function controls NNLO contributions in resummed 
calculations at full next-to-next-to-leading logarithmic accuracy. The same coefficient 
function is used in applications of the subtraction method to perform fully-exclusive 
perturbative calculations up to NNLO. 
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The transverse-momentum (qx) distribution of systems with high invariant mass M (such 
as Drell-Yan lepton pairs, photon pairs, vector bosons, Higgs bosons, and so forth) produced 
in hadron collisions is computable by using perturbative QCD. However, in the small-g^ region 
(roughly, in the region where qr <C M) the convergence of the fixed-order perturbative expansion 
in powers of the QCD coupling as is spoiled by the presence of large logarithmic contributions 
of the type ln n (M 2 /gf,). The predictivity of perturbative QCD can be recovered through the 
summation of these logarithmically-enhanced contributions to all order in a§ pQ. 

The structure of the resummed calculation can be organized in a process-independent form 
[21 El IU [5] , in which the logarithmic contributions are controlled by a set of perturbative functions, 
usually denoted as A(as), B(as), C(as) and H(cns) ( see 5 e -g- ; Eqs. (jSJ) and (fTBl) and related com- 
ments). These functions and, hence, their perturbative coefficients (e.g. the coefficient of the 
n-th order contribution A^a§ to A(ag)), have no explicit dependence on the ratio qx/M. The 
perturbative coefficients, once they are known, can be inserted in process-independent resumma- 
tion formulae that systematically resum, in explicit form, the classes of leading, next-to-leading, 
next-to-next-to-leading (and so forth) logarithmic contributions to the transverse-momentum dis- 
tribution. In this respect, the transverse-momentum resummation program has formal analogies^ 
with the study of logarithmic scaling violations (of ultraviolet or collinear origin), where the re- 
summation of logarithmic terms is traded for the calculation of perturbative functions, such as 
short-distance coefficient functions and anomalous dimensions. 

Most of the qT resummation coefficients are known, since some time 0, [71 [HI [9] , up to the second 
order in ag. The third-order coefficient A^ has been obtained in Ref. [10]. In recent years, we 
have been working on a research project aimed at the completion of the qT resummation program 
at the second perturbative order. This requires the calculation of the second-order coefficient 
function^ 'H^(z) (see Eq. (fl9]) ). which includes a pro cess- dependent part. The computation of 
the coefficients has been explicitly carried out for two benchmark processes, namely, Higgs 
boson production and the Drell-Yan process, and the corresponding results have been obtained 
and used in Refs. [IT] and [12], respectively. In the case of the Drell-Yan process, the result has 
also been applied [J3] to the qx spectrum of the Z boson, by explicitly performing transverse- 
momentum resummation at the next-to-next-to-leading logarithmic (NNLL) accuracy. In this 
paper, we consider Higgs boson production and we document the result for [11] in explicit 
analytic form. We also illustrate the method that we have used to perform the calculation. 

Considering the perturbative contributions that are logarithmically-enhanced and, thus, sin- 
gular in the limit q<r — > 0, the calculation of completes the knowledge of the qT distribution 
at full next-to-next-to-leading order (NNLO) accuracy. This fact has implications in the context 
of both resummed and fixed-order calculations, as we briefly discuss below. 

Some recent resummed calculations of the qx spectrum of the Standard Model (SM) Higgs 
boson at Tevatron and LHC energies are presented in Refs. [H]-[2I]- The inclusion of in 
calculations that use the other qj- resummation coefficients up to NNLL order gives theoretical 
predictions that embody the exact NNLO calculation in the small-g;r region. The NNLL resummed 
calculations can then be properly matched (by using, for instance, qx resummation as in the 

^These analogies may hide important physical, conceptual and technical differences, which are discussed in the 
literature on transverse-momentum resummation. 

•tin this introductory part we are using a shorthand notation, since the symbol H^ 1 actually refers to a set of 
several coefficient functions. 
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formalism of Ref. [T9]) with the customary fixed-order calculation at large q?, in such a way that 
the integration over q^ of the qx distribution exactly returns the NNLO value of the total cross 
section. Indeed, a rough approximation of 'HS 2 \ such as to reproduce the NNLO value of the 
total cross section with good numerical accuracy, was constructed and used in Ref. [19]. The 
approximation of Ref. [19] represents a very crude estimate of the function W?\z)\ nonetheless, 
that approximation quantitatively works very well (especially at LHC energies) pJJJ [20] over a wide 
range of Higgs boson masses. An updated version of the code HqT [19J, which also implements the 
exact coefficient T-L^ 2 \ will be available shortly. 

In Ref. [UJ, we presented a practical formalism to perform NNLO calculations at the fully- 
exclusive level for a specific class of processes, namely, the production of colourless high-mass 
systems in hadron collisions. The formalism exploits the subtraction method to cancel the unphys- 
ical infrared divergences that separately occur in the real and virtual radiative corrections. The 
explicit construction of the subtraction counterterms [TU [20] is based on the process-independent 
structure of transverse-momentum resummation formulae and on their expansion up to NNLO in 
QCD perturbation theory. The formalism thus requires the complete knowledge of the qr resum- 
mation coefficients up to 0(a\). Although the results of the present paper were not explicitly 
illustrated in Ref. [TT], they were taken into account in the NNLO computations presented therein. 
In particular, the explicit application to Higgs boson production (which was implemented in the 
Monte Carlo code HNNLO) considered in Refs. [TTj [22] is based on and implements the analytic 
results for the coefficient function that are documented in the present paper. 

The paper is organized as follows. We first introduce our notation and describe the small- 
er behaviour of the Higgs boson cross section up to NNLO. Then we briefly review transverse- 
momentum resummation for Higgs boson production and the corresponding all-order resummation 
formula recently derived in Ref. [5]. The new resummation formula differs from its naive version 
that is used in the literature: the differences start at C?(a|), which is relevant for the purposes 
of the present paper. Finally, we present the analytic results of our NNLO calculation of the gr 
distribution. The results are expressed directly in terms of the coefficient function 'H^(z) and 
related resummation coefficients. We conclude the paper by describing the method that we have 
used to perform the NNLO calculation. 

We briefly introduce the theoretical framework and our notation. We consider the production 
of the SM Higgs boson H, through the gluon fusion mechanism gg — > H, in hadron-hadron 
collisions. The effective coupling ggH is produced by heavy-quark loops, and the top quark 
gives the dominant contribution. We treat the coupling ggH in the framework of the large-m top 
approximation [231 [2U and we consider a single heavy quark, the top quark with mass m top , 
and np (jif = 5) massless- quark flavours. We use the narrow width approximation and we treat 
the Higgs boson as an on-shell particle with mass M. The QCD expression of the Higgs boson 
transverse-momentum cross section is 

do~ do 

-TT(q T ,M,s) = J2 dz i / dz 2 f a/h Az u M 2 )f b/h Az 2 ,M 2 )-^(q T ,M 1 s = z 1 z 2 s-,a s (M 2 )), 

«<?T ab JO JO "Qt 

(1) 

where f a /hi(%, Hf) ( a = d) are the parton densities of the colliding hadrons {hi and h 2 ) at the 

factorization scale fip, and da ab /dq T are the partonic cross sections. The centre-of-mass energy 
of the two colliding hadrons is denoted by s, and s is the partonic centre-of-mass energy. We 
use parton densities as defined in the MS factorization scheme, and as(fi 2 R ) is the QCD running 
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coupling at the renormalization scale in the MS renormalization scheme. In Eq. ([T]) and 
throughout the paper, the arbitrary factorization and renormalization scales, fip and [Mr, are set 
to be equal to the Higgs boson mass M. 

The partonic cross sections da a bjdq\ are computable in QCD perturbation theory as power 
series expansions in as{M 2 ). We are interested in the perturbative contributions that are large in 
the small-gr region (q? <C M) and, eventually, singular in the limit qx — > 0. To explicitly present 
the perturbative structure of these enhanced terms at small gr, we integrate the gr distribution 
over the region < gr < Qo, and we introduce the cumulative partonic cross section 

[ Q °dq 2 T ^(q T ,M,s = M 2 /z;a s (M 2 )) = za^(a s (M 2 )) R ab (z, M/Q ;a s (M 2 )) , (2) 
Jo ^Qt 

where the overall normalization of the function R ab is defined with respect to cr^, which is the 
Born level cross section for the partonic subprocess gg — > H. Using the large-mt op approximation, 
the explicit expression of is (23] 

{as) = 5iwi • (3) 

where Gp is the Fermi constant. The partonic function R has the following perturbative expansion 

oo 

R ab (z,M/Q 0] a s ) = 5 ga 5 gb 5(l - z) + ^ (^)™ R^(z,M/Q ) . (4) 

n=l 

The next-to-leading order (NLO) and NNLO contributions to the cumulative cross section in 
Eq. ([2]) are determined by the functions Iv- 1 ' and R^ 2 \ respectively. The small-gr region of the 
cross section da a bldq\ is probed by performing the limit Qo <C M in Eq. ([2]). In this limit, the 
NLO and NNLO functions R^ and R^ have the following behaviour: 

R$(z, M/Q ) = l 2 B$*\z) + k + RT\z) + O(Q 2 /M 2 ) , (5) 



R${*> M /Qo) = it R%%) + i 3 o RT\z) + il RT(z) + k R% 1] (z) + RT(z) + O(Q 2 /M 2 ) , 

(6) 

where l = lia(M 2 /Ql). In Eqs. ([5]) and ([6]), the powers of the large logarithm Z are produced 
by the singular (though, integrable) behaviour of da a bjdq\ at small values of g^. The coefficients 
rO-'' 171 ) (with m < 2) and R^ 2]m ^ (with m < 4) of the large logarithms are independent of Qo', 
these coefficients depend on the partonic centre-of-mass energy s and, more precisely, they are 
functions of the energy fraction z = M 2 /s. 

In this paper we present the computation of the cumulative cross section in Eq. (T5]) up to 
NNLO. The partonic calculation is performed in analytic form by neglecting terms of 0(Ql/M 2 ) 
in the limit Qo *C M. Therefore, we determine the coefficient functions R( n]m \z) in Eqs. ([5]) and 
(Q. Before presenting the results, we illustrate how these functions are related to the perturbative 
coefficients of the transverse- momentum resummation formula for Higgs boson production [3, 5J. 
This relation, which allows us to extract the g^ resummation coefficients up to O(ag), also shows 
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that the knowledge of Eq. (E]) is sufficient to fully determine the NNLO rapidity distribution of 
the Higgs boson in the small-g^ region. 

The partonic cross section da a b/dq T in Eq. (pQ) can be decomposed in the form da a b = da^ b mg ^ + 
d&ab- The singular component, da^ b mg \ contains all the contributions that are enhanced at 
small qx- These contributions are proportional to 5(q T ) or to large logarithm^) of the type 
4r \n m (M 2 / q T ) . On the contrary, the remaining component, da* b , of the partonic cross sec- 
tion is regular order-by-order in a s as qx — > 0. To be precise, the integration of da^ b eg) /dq T over 
the range < qx < Qo leads to a finite result that, at each fixed order in as, vanishes in the 
limit Qo — > 0. Therefore, da^ b eg ^ only contributes to the terms of 0(Qq/M 2 ) on the right-hand 
side of Eqs. ([5]) and ([6]). The decomposition of the partonic cross sections can be inserted in the 
right-hand side of Eq. (TjQ), thus leading to the corresponding decomposition of the hadronic cross 
section da /dq T , namely, da = da^ 111 ^ + da^ reg \ 

We consider the singular component of the Higgs boson qx cross section, and we recall its 
all-order perturbative structure. We directly refer to the hadronic cross section (rather than the 
partonic cross sections), since its structure can be presented by using a more compact notation. 
Moreover, to illustrate the general kinematics of transverse-momentum resummation, we consider 
the qx cross section da/dydq T at fixed value of the rapidity y of the Higgs boson (the rapidity 
is defined in the centre-of-mass frame of the two colliding hadrons). The transverse- momentum 
resummation formula for the singular component of the Higgs boson cross section is [U [5] 



da {sing ) , , M 

dy dq, 



(y, q T , M, s) — ag } (as(M 2 )) / db - J Q (bq T ) S g {M, b) 



X 

11, «2 



E f— f— [ HFC ^] 99 -a ia2 fa l/h Axi/zubl/b')f a2/h2 ( X2 /z 2 X/b 2 ) , (7) 



where the kinematical variables x« {i = 1, 2) are X\ = e +y M/ y/s and x 2 = e~ y M/ y^s. The integra- 
tion variable b is the impact parameter, Jo(bqx) is the Oth-order Bessel function, and bo = 2e~ lE 
(jE = 0.5772... is the Euler number) is a numerical coefficient. The symbol \H F CiC2 1 \ gg . aia2 
denotes (see Eq. (44) in Ref. jS]) the following function of the longitudinal-momentum fractions 
Zi and z 2 : 



[H H CiC 2 ] gg . aia2 = H g H (a s (M 2 )) [ C 9ai { Zl -a s (bl/b 2 )) C ga , 2 (z 2] a s (b 2 /b 2 )) 

+ G gai (z 1 ;Mbl/b 2 )) G ga2 (z 2 ;a s (b 2 /b 2 ))} , (8) 



where if^(as), C ga (z; as) and G ga (z; as) are perturbative functions of as (see Eqs. (!9l- (jTT!) ). The 
other perturbative ingredients of Eq. ([7]) are the function S g (M,b), which is the Sudakov form 
factor of the gluon (see the comments below), and a^\as(M 2 )), which is the Born level cross 
section in Eq. ([3]). 

The structure of Eq. ([7]) is well known in the literature on resummed calculations for the qx 
spectrum of the Higgs boson. However, the functional form of Eq. ([8]) is new [5]. The customary 



§To be precise, the logarithms are combined with corresponding 'contact' terms, which are proportional to S(q^). 
These combinations define regularized (integrable) 'plus distributions' ^- ln TO (M /gf.) with respect to q^. The 
cumulative cross section in Eq. ^) is insensitive to the precise mathematical definition of these 'plus distributions'. 
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'naive' version (i.e. the version extrapolated from the transverse-momentum resummation formula 
for the Drell-Yan process) of Eq. (JSJ) includes only the perturbative functions and C ga . The 
presence of an additional term, due to the function G ga , has been pointed out in Ref. [5]. Note 
that the function G ga (z; as) is of 0(ats) ( see Eq. (j2j)) and, therefore, it leads to a contribution of 
0(ais) in the right-hand side of Eq. (JBj). This fact implies that the presence of G ga (z; as) cannot 
be detected through a NLO calculation of the qr spectrum of the Higgs boson. This fact also 
implies that our NNLO analytic calculation of the cumulative cross section in Eq. (j2J) gives an 
explicit check of the presence of G ga (z;as) and of its precise form at 0(as) (see Eq. (1281) and 
related comments). 

The gluon form factor S g (M,b) of Eq. (J7j) is a process- independent quantity [3J, El SJ- Its 
functional dependence on M and b is controlled by two perturbative functions, which are usually 
denoted as A g (as) and B g (as) (see, e.g., Eqs. (10)-(12) in Ref. [5]). Their corresponding n-th 
order perturbative coefficients are A g n ^ and B g n \ The coefficients A g , B g l \ A g |8J and B g [9] 
are known: their knowledge fully determines the perturbative expression of S g (M, b) up to 0(a§). 

The quantity [if F CiC 2 ] in Eq. (J7|) depends on the three perturbative functions Hg ,C ga and 
G ga . By inspection of the right-hand side of Eq. (jSJ), we notice that the scale of as is not set 
to a unique value. We have as(M 2 ) in the case of the function Hj^(as), and asib^/b 2 ) in the 
case of the functions C ga (as) and G ga (as)- The presence of these two different arguments of as is 
related to the physical origin [H [5] of the corresponding perturbative functions. Roughly speaking, 

\as(M 2 )) embodies contributions due to the /iard-momentum regiorflof the virtual corrections 
to the lowest-order subprocess gg —> H. The functions C ga and G ga instead refer to the inclusive 
subprocess ga — > H + X: roughly speaking, C ga (as(bl/b 2 )) and G ga (as(bl/b 2 )) originate from 
the kinematical region where the momenta of the partons in the final-state system X are (almost) 
collinear to the momentum of the initial-state parton a. Owing to this physical picture, the 
quantity [if^CiC^] can be regarded as a hard- collinear partonic function. Note that the function 

(as) is process dependent, since it is directly related to the production mechanism of the SM 
Higgs boson. On the contrary, the partonic functions C ga and G ga are process independent, as a 
consequence of the universality features of QCD collinear radiation. 

We recall that the functions (as),C ga (as),G ga (as) and the perturbative function B g (as) 
of the gluon form factor are not separately computable in an unambiguous way. Indeed, these four 
functions are correlated (constrained) by a renormalization-group symmetry [I] that is related 
to the 6-space factorization structure of Eqs. (J7j) and (jSj). The unambiguous definition of these 
four functions thus requires the specification^ of a resummation scheme [1] . Note, however, that 
considering the perturbative expansion! of Eq. (J7|) (i.e., the perturbative expansion of the singular 
component of the qr cross section), the resummation-scheme dependence exactly cancels order- 
by-order in as- 

The perturbative expansion of the three functions on the right-hand side of Eq. (JSJ) is defined 

^This is the region where the size of the momenta of the virtual loops is of the order of M. 

'■The reader who is not interested in issues related to the specification of a resummation scheme can simply 
assume that Hg(ots) = 1 throughout this paper. 

§The resummation-scheme dependence also cancels by consistently expanding Eq. (T7J) in terms of classes of 
resummed (leading, next-to- leading and so forth) logarithmic contributions [19] . 
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as follows: 

oo 

G 9a (z;a s ) = ^GU(z)+J2{^yG<?Kz) , (9) 

n=2 



C ga ( Z] as) = 5 ga 5(l-z) + J2(^yci n Xz) , (10) 

n=l 

oo 

ff>s) = i + E(v)"^ < " ) - < u > 



7T 

n=l 



Since the partonic functions G ga and C ga are process independent, they fulfil the following rela- 
tions: 

G gqf (z;a s ) = G g g fl (z;a s ) = G gq (z;a s ) , C gqf (z; a s ) = C g g f ,(z; a s ) = C gq (z; a s ) , (12) 

which are a consequence of charge conjugation invariance and flavour symmetry of QCD. The 
dependence of G ga (C ga ) on the parton label a is thus fully specified by G gg and G gq {C gg 
and C gq ). The first-order coefficient functions G gg (z) and G gq (z) (they are independent of the 
resummation scheme) are known [5]: 

Gfl{z) = C A , GJJC*) = C F . (13) 

The first-order coefficient function Cgq(z) is also independent on the resummation scheme; its 
expression is [26] 

Cfl{z)= l -C F z . (14) 

Using the large-m top approximation, the first-order coefficients H^ 1 ^ and Cgg (z) fulfil the follow- 
ing relation [H] : 

+ ^ «(1 " *) = (5 + 7F2) ^' 3Cf *(1 - z) • (15) 

The separate determination of C g g\z) and H^ 1 ^ requires the specification of a resummation 
scheme. For instance, considering the resummation scheme in which the SM Higgs boson coefficient 
vanishes, the right-hand side of Eq. f|T5|) gives the value of Cgg(z) [26], and the corresponding 
value of the gluon form factor coefficient B g is explicitly reported in Eq. (128) of the second paper 
in Ref. [9]. The computation of the second-order coefficients Cgq, C g 2 g and H^ 2 ^ is the aim of the 
calculation described in this paper. 

For later purposes, we also define the following hard-collinear coefficient function: 

Kfg*-ab( z > a s) = Hg (as) / dz\ \ dz 2 6(z-z 1 z 2 ) C ga (z 1 ;as)C gb (z2;as)+G ga (z 1 ;as)G gb (z 2 ;as) 

Jo Jo 1 

(16) 

which is directly related to the coefficient function in Eq. (jSJ). There are only two differences 
between Eqs. (jSJ) and ffTB"]) . The first difference is due to the fact that the function % H depends 
on the energy fraction z, since the right-hand side of Eq. ffTB"]) involves a convolution integral 
over the momentum fractions Z\ and z 2 . This convolution kinematically arises by considering the 
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integration of Eq. (J7J) over the rapidity y of the Higgs boson. The second difference regards the 
scale of as: in the functions H H (as), C(as) and G(a$) on the right-hand side of Eq. (Tl6l) . the 
argument of as is set to the same value (this common scale is not explicitly denoted in Eq. (fT6j) ) . 
Owing to this feature, the process-dependent function T-L^ g ^ ab is unambiguously defined (i.e., it is 
independent of the specification of the resummation scheme) [I]. The perturbative expansion of 
the function TL directly follows from Eqs. (I9l)-(|TT1). We have: 

oo 

K^(z-,* S ) = 5 ga 5 gb 5(l-z)+J2(^ Ytf^Jz) , (17) 

n=l 

where the first-order and second-order contributions are 

= S ta 6 gb S(l - z) flf « + 5 ga cf h \z) + 5 gh Cfl{z) , (18) 

= 8 aa 6 Bb 5(l - z) Hf^ + 8 ga Cf b \z) + 5 gb Cfl{z) + tff « (s ga C { g l b \z) + 5 gb C^{z)) 
+ ® C S) (*) + ® G$) (*) . (19) 

In Eq. (TX9T) and in the following, the symbol £g> denotes the convolution integral (i.e., we define 
(g <g> h)(z) = J dzi Jq 1 dz 2 S(z - z 1 z 2 ) g{zx) h(z 2 )). 

After our illustration of the all-order resummation formula in Eq. (JTJ), we can return to its 
relation with the perturbative expression of the cumulative partonic cross section in Eq. (j2J). Using 
the Altarelli-Parisi evolution equations, the parton densities f a /h(x,bl/b 2 ) on the right-hand side 
of Eq. (j7|) can be expressed in terms of the corresponding parton densities f a / h (x, M 2 ) at the 
factorization (evolution) scale fip = M. Having done that, all the remaining factors in Eq. ([Tj) 
are the partonic contributions that determine the small-gy singular component of the Higgs boson 
partonic cross section da ab /dq^ in Eqs. (JTJ) and (J2]). At fixed values of the impact parameter b, all 
these partonic contributions can be expanded in powers of as(M 2 ), thus leading to perturbative 
coefficients that depend on powers of ln(6 2 M 2 ). The dependence on ln(6 2 M 2 ) is produced by the 
gluon form factor S g (M, b), by the Altarelli-Parisi evolution equations and by the QCD coupling^] 
«s(&o/^ 2 )- The powers of ln(6 2 M 2 ) can then be transformed into logarithms, ln(M 2 /q^), in qT- 
space by explicitly performing the Bessel transformation (i.e. the integration over b) in Eq. ([7]). 
This procedure, which involves manipulations that are standard in the context of transverse- 
momentum resummation (technical details can be found, for instance, in Ref. [IS]), yields the 
explicit perturbative expression of the singular component of da ab j dydq\ . The integration of this 
expression over y (y is the rapidity of the Higgs boson in the centre-of-mass frame of the two 
colliding partons) and qx finally gives the cumulative partonic cross section of Eqs. (|2J) and (jlj) 
in the limit Q < M (i.e., modulo the contributions of 0(Ql/M 2 ) in Eqs. (ED) and ©). The 
perturbative functions R^\z,M/Q ) of Eq. (jl]) have a dependence on l = \n(M 2 /Ql) that is 
explicitly determined by the resummation formula (j7j), whereas the dependence on z is given in 
terms of the qx resummation coefficients (those of the gluon form factor and in Eqs. (J9"j)- (|11|) ). 

The NLO and NNLO functions R^ b and R^ b have the following expressions: 

R { 2 (*, M/Q ) = ll £f 9 S(z) + l Sf^(z) + H^Jz) + 0{Q 2 JM 2 ) , (20) 

"'"The coupling as{b^/b 2 ) can be expressed in terms of a^{M 2 ) and ln(6 2 M 2 / '6q) by using the renormalization 
group equation for the perturbative /Lt 2 -evolution of the running coupling a^{^ 2 )- 
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+ (<2*M - <3 £f 9 22W) + 0(QS/M 2 ) , (21) 

which are consistent with the behaviour in Eqs. (jSJ) and dHj). In Eqs. (1201 and (l2Tj) we use the 
same notation as in Ref. [13]. The coefficient functions ^fg^ab(z) depend on the resummation 
coefficients: the explicit expressions are given in Eqs. (63),(64),(66)-(69) of Ref. [19] (we have 
to set Hr = /jlf = Q = M, where /j,r, (jlf and Q are the auxiliary scales of Ref. [19J) and 
are not reported here. The coefficients H^g^ ab and 1-L^^_ ah are exactly those in Eqs. ( TT8l) and 
( TT9i) . The first-order terms and depend on the gluon form factor S g (M,b). The 

second-order terms depend on %^_ a6 and on the gluon form factor S g (M, b) up to C(ag). 

The numerical coefficient £3 ~ 1.202 ... is the Riemann ^-function) on the right-hand side of 
Eq. (12T1) originates from the Bessel transformations (see, e.g., Eqs. (B.18) and (B.30) in Appendix B 
of Ref. US]). 

The relations (|20|) and (12 ip can be exploited in two different ways. From the knowledge of the 
perturbative coefficients of the resummation formulae (J7J) and OH]), we can compute Y, H ( n ' m ) and 
ft H < n > and then, we can obtain a perturbative prediction for the cumulative partonic cross section 
up to NNLO. Alternatively, from the explicit NNLO perturbative computation of the cumulative 
partonic cross section, we can extract £^ m ) and U H{n) and then, we can determine the qx 
resummation coefficients up to 0(a\). 
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Our NNLO computation of the cumulative partonic cross section is described in the final 
part of this paper. We obtain the following results. The explicit result of the NLO function 

(z), Yi^'lliz) and T-L^„^ nh (z), as predicted by the q F 



d { 2(z) confirms the expressions of E^j2)(s), ^ g tab( z ) and Uf^y 
resummation coefficients at 0(as). At NNLO, the present knowledge [H [9] of the qx resummation 
coefficients at 0(a$) predicts the expressions of the terms S^^(z), with m = 1,2,3,4. Our 
result for the NNLO function R^(z) confirms this prediction, and it allows us to extract the 
explicit expression of the second-order coefficient function %^g^ a 



-ab\ 



We obtain: 



2(1 



(2 + z)< 
4z 



In z 



(22) 



= Cl\ i (2 - z) In 3 z - ^ (3.- + 4) In 2 z + ^ (z - 3) In z 



+ M^H ln3(i - z)+ sH- 6)ln2(i - z) 

+ Af + ?_ 2 V (1 _ z) + 5_H z 

V 8 Z J 8 16 

+ C F n F LL (1 + (1 - z) 2 ) ln 2 (l - z) + 1 (z + I - 5) ln(l - 



14 14 13 

1 1 z 

27 27z 108 



+ CfCa i - ii±4±^L i3 m + fLi-^i Li3W 



2z V 1 + z I \2 2z 4 



3 / , > 2 \ T / \ f 11 z z 2 ( 1 3 3z\, \ 
_(l + ( 1 + ^)L,3(-,) + (2---- + 3 + (-- + - + T jln,jLi 2W 



,'z (l + (l + z) 2 ), .... 
(4-- ^^hi(z))Li 2 (-z) 



+ (1 + ^ ln 3 (l + «) + ^ ( (1 + (1 + zf) (3 In 2 z - tt 2 ) - 6z 2 In 0) ln(l + z) 

~ (1 + ( 2 1 4 7 )2)) m 3 (l -z) + ± (6(1 + (1 - ^) 2 ) In, - 5z 2 - 22(1 - z)) ln 2 (l - z) 

+ — (-152 + 152z - 43z 2 + 6(-22 + 24z - 9z 2 + 4z 3 ) In 2 + 9(1 + (1 - z) 2 ) In 2 z) ln(l - 



107 1 z 11,,, 

H z 2 In 2 



- — f 1 + In 3 z + — (36 + 9z + 8z 2 ) In 2 z + ( 

12 V 2/ 48 V ; V 24 * 12 9 

1/, 503 11 ,\ 1007 tt 2 5, /tt 2 n> 133\ 2 /38 tt 2 

+ - 4C 3 + — vr 2 + ( 3 + z[ — + 2( 3 + z 2 

zl 54 36 J 108 3 2 s3 I 3 108/ V 27 18 



(23) 
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1 A V 288 8 72 144 18 S 7 V 24 8 4 

9 ^2 5^ \„ ^ /287 5 2 4 ' 

+ -Cl C A C F -C A n F \ 1 7T 2 + -C 3 

4 F 96 12 V 144 36 9 s 

., , (1 + z + z 2 ) 2 / f z \ T , A 2 - 17^ - 22z 2 - 10z 3 - 12z 4 > 
+ C 2 ^ .., ■ / 2Li 3 — — -Li 3 (-z) + - ■ . C: 



z(l + z) V \ l + z J J 2z(l + z 



3 



5 - z + 5z 2 + z 3 - 5z 4 + z 5 /T , . , N . . ln(» 3 - z + 3z 2 + z 3 - 3z 4 + z 5 
z(l-z)(l + z) (L13(Z) " Cs) + Ll2( " } — zJTTz) 

+ ^^l^ 2 (>{z)U 2 {-z) - i ln 3 (l + z) + C 2 ln(l + z] 

l-z, 9 . r . 1 , . , 1 ln 3 (z) (1 + z - z 2 ) 2 

, 2 , . /(I-2 + 2 2 ) 2 , „ , (1 + , + r 2 ) 2 , „ . 25-ll, + 44, 2 
+ h « ( 2,(1-,) - Z) - 2,(1 + ,) ln(1 + Z) + 24 

, , V(l + , + 2 2 ) 2 , 2 ,, . (l-, + 2 2 ) 2 , 2 ,, . 

72 + 773z + 149z 2 + 536z 3 \ 517 449 380z 835z 2 ' 



72z J 27 27z 27 54 

r , „, I i±£ i n 2 (z) + 1(13 + lOz) ln(z) - - ln(l - z) - - + — + -z - — z 2 
1 12 v ; 36 V ; v ; 12 v ; 54 108z 54 108 

+ C F n F \^ ln 3 (z) + ~(3 + z) ln 2 (z) + 2(1 + z) ln(z) - ±^(1 - 23z + z 2 ) 

(24) 

where L t = ln(M 2 /m 2 op ) (m top is the pole mass of the top quark) and Li fc (z) (k = 2,3) are the 
usual polylogarithm functions, 

r z dt r l dt 

Li 2 (z) = - / T ln(l -t) , Li 3 (z) = / a - ln(t) ln(l - zt) . (25) 



We comment on the Higgs boson results in Eqs. (j22p - (|24p and on the ensuing determination 
of the second-order coefficients Cg 2 q,Cg 2 g and H^ 2 ^ in Eqs. (fT0|) and fTTTjl . 

Considering the dependence on the parton indeces a and b, the hard-collinear function ?/^ < _ afe 
is 'flavour blind', namely, it fulfils the relations 

-T/# _ -T/ff _ -T/H _ -T/W — 'JjH ('261 

99<-QfQfi 99^-QfQfi 99^-QfQf 99<-Qf9f ~ 99^11 ' V > 
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njH _ njH _ njH _ niH _ njH (271 

99<r-1f9 99^9Qf 99^9Q f ' 99<~5/'9 ~ 99^91 ' N > 

These symmetry relations follows from Eq. (1121) and from the convolution integral in Eq. (TT6j) 
(the convolution integral implies that T-tg gi _ ab is symmetric with respect to the exchange a -H- b). 

Therefore, the parton matrix Tifg^j, is completely specified by the three entries in Eqs. (122|) -( |24|) : 

the quark-quark function 'H^J^qq, the gluon-quark function H^g^gq and the gluon-gluon function 

n gg^gg- 

Using Eq. ( !T9|) . in the quark-quark channel we have 

nf g % q (z) = (C7W ® <7«) (z) + (G« ® GW) (z) . (28) 

We see that the second-order coefficient function %^i- 3g (z) is fully determined by the resumma- 
tion coefficients at 0{a$). Using the values of Gg\ and Cgq in Eqs. ( TT3|) and (fT4|) . the expression 
on the right-hand side of Eq. (1281 is in complete agreement with the result in Eq. (!22l . Therefore, 



our explicit computation of the NNLO partonic function R qg represents a consistency check of the 
factorization formula (jBJ) and of the value of Ggq, which were derived in the process-independent 
study of Ref. [5]. 

Considering the gluon-quark channel, Eq. (I19p can be recast in the following form: 

Cfliz) + I HfV C«(z) = H? g % q (z) - \ (HfMgg ® C$) (z) - (G« ® G«) (z) , (29) 

where we have used Wfg^gg(z) = Hg 5(1 — z) + 2G£g(z) (see Eq. ( {TBI ). The relation ([291 can 
be used to determine Cfq\z) from the knowledge of Hgg^gq and of the gr resummation coefficients 
at 0(ag) (in particular, the values of Ggg and Ggq affect the determination of Cgq). Inserting 
the first-order results of Eqs. (fl3"l) - (fT51) in Eq. ([2"9l . we explicitly have: 



1 + z) In z + 2(1 -z)- l±*L z 2 



+ \ < W ^ ^ = K ( ^) + Cl\z + C F C A \ 

"(30) 

where Ugg^ gq is given in Eq. fl23]). Note that the right-hand side of Eq. (J29]) (or Eq. fl3Qj)) 
is resummation-scheme independent. Analogously to Eq. (|T5|) . the dependence of Cg ? on the 
resummation scheme is thus parametrized by the first-order coefficient H g on the left-hand side 
of Eq. (ED]). 

The process-independent coefficient function Cg 2 g \z) is obtained analogously to Cg 2 q \z). Con- 
sidering the gluon-gluon channel, Eq. ( Tl9|) gives: 



2 <?<?(*) + *(!"*) 



g a V f / 



"gg<^gg^ z ) 



, 1 o-if(i) 

where the right-hand side is expressed in terms of resummation-scheme independent functions. 
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Inserting Eqs. ( 1T3"|) -( 1T5|) in Eq. (I3"T|) . we explicitly obtain: 



H H <$ - - (H H ^) 2 + ( 5 + k 2 )C a - 3C F H(1) 

9 4 V 9 / ~ j 9 



2 C?^) + 



+ 



(5 + tt 2 )C a - 3C F 



2 1 



+ C 2 A -[(l + z)\nz + 2(l- z)] , (32) 



z 



where T-L gg \J gg is given in Eq. ( 1241) . We observe that Cgg(z) includes a resummation-scheme 
dependent part that is simply proportional to 5(1 — z). This part depends on and H^ 2 \ 

We also recall [4] that the resummation-scheme invariance relates Cg g , if <f^ 2 ^ and the third-order 

(3) 

coefficient B g of the gluon form factor. 

We conclude the paper by briefly describing the method that we have used to perform the 
NNLO analytic computation of the Higgs boson cross section in Eq. (j2J). The NNLO partonic cal- 
culation has to be carried out by using dimensional regularization to evaluate the QCD scattering 
amplitudes and their integration over the partonic phase space. In the framework of the large-m top 
approximation, the relevant partonic subprocesses are: the gluon fusion subprocess gg — >■ H up 
to the two-loop level, the single-emission subprocesses ab — > H + c up to the one-loop level, and 
the double-emission subprocesses ab — > H + c\ + C2 at the tree level. The corresponding scattering 
amplitudes are known and have been used in the analytic calculations of two relevant Higgs boson 
observables: the NNLO total cross section o^f* [27l[28| [29] and the NLO differential cross section 
da a b/dydq^ [301 E] at large qr- To perform our NNLO calculation, we take advantage of these 
available results: both observables are indeed computed up to relative order «g with respect to 
the Born level cross section (as)- We rewrite the integration in Eq. ([2]) as follows: 

/ d lT -TT{QT,M,s;a s ) = [ dq? ^^{q T ,M,s;as) - [ dq\ ^^(q T , M, s\ a s ) 
Jo dq^ J Q dq£ J Q 2 dq^ 

= d-vJ rt )(JVf,s;a s ) - / dq% \ dy f (y,q T ,M,s;a s ) . (33) 
JqI J-oo dy dq z T 

The cumulative partonic cross section over the range < qr < Qo is thus obtained by subtraction 
according to Eq. ( |33l) : we start from the total cross section!] <x^° t ' ) and we subtract the contribution 
due to the qr cross section in the 'large '-gr region where qx > Qo (in the context of Eq. ( 1331) . 
'large' values of qx generically means 'non- vanishing' values of qx). The differential cross section 
dcrabl 'dy dq\ in the integrand on the right-hand side of Eq. (133!) is presented in Ref. [30] in complete 
analytic form: we use this form and we explicitly carried out the integrations over y and qx- Since 
Qt > Qoi these integrations can directly be performed in four space-time dimensions, with no 
further use of dimensional regularization. At NLO, the cumulative partonic cross section can be 
computed in explicit analytic form0 for arbitrary values of Qo- At the NNLO, we limit ourselves 
to analytically computing the cumulative cross section in the limit Qo ^ M, thus neglecting the 
terms of 0(Ql/M 2 ) on the right-hand side of Eqs. (jSJ) or (|2"T|) . 



§ We actually use the expressions of Ref. [29] . which are given for general colour factors, Cf — (AT? — l)/(2iV c ) 
andCU = N c , of SU(N C ). 

^The NLO analytic result was already reported in Sect. 4 of Ref. [32]. Indeed, at the NLO, the cumulative 
cross section of Eq. ([2]) exactly coincides with the jet- vetoed cross section cr vcto of Ref. [32] , provided we identify 
Qo = P~r to ', where p™ to is the jet veto parameter. 
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In this paper we have considered the production of the SM Higgs boson at hadron colliders. 
We have presented the NNLO analytic calculation of the cross section at small values of g-r 
(see Eqs. (j2J) and (jSJ)). The NNLO result is compared (see Eq. ff2~Tj) ) with the predictions of 
transverse- momentum resummation (see Eqs. ([7]) and (jEJ)). The comparison gives a crosscheck 
of the factorization formula (jHJ) (see also Eq. ( l28i) ) and allows us to determine the previously 
unknown resummation coefficients at 0{a§). These are the coefficient functions 'Hfg^l a b(z) (see 

Eqs. (j22l-(j24l) and the related coefficients Cf] and Cf] (see Eqs. (j30l) and (j32))). which control the 
dependence on the rapidity of the Higgs boson. These coefficients can be implemented in resummed 
calculations of the inclusive qx distribution at full NNLL accuracy. Using the method of Ref. |11| . 
the same coefficients are necessary to perform the fully-exclusive perturbative calculation up to 
NNLO. 
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